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We use the Bethe Ansatz solution for the one dimensional Hubbard model with open 
^ . boundary conditions and applied boundary fields to study the spectrum of bound states 

C^ ' at the boundary. Depending on the strength of the boundary potentials one finds that the 

true ground state contains a single charge or, for boundary potentials comparable to the 
Hubbard interaction, a pair of electrons in a bound state. If these are left unoccupied one 
finds holon and spinon bound states. We compute the finite size corrections to the low 
lying energies in this system and use the predictions of boundary conformal field theory 
to study the exponents related to the orthogonality catastrophe. 
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1 Introduction 

The recent advances in the understanding of boundary effects in low dimensional quantum systems 
due to the predictions of boundary conformal field theory [||, |2|, y] and the formulation of Bethe 
Ansatz soluble models on open lattices with potentials applied on the boundary sites 0, |5[ g have 
opened new possibilities to study the effects of correlations and quantum fluctuations on long standing 
problems such as the orthogonality catastrophe P, |9[ and edge singularities in optical absorption 



experiments [|10|, |11|, |12| . 

The effect of electronic correlations on the bulk critical behaviour of (1 + l)-dimensional quantum 
systems has been studied successfully in the Tomonaga-Luttinger model which then can be handled 



using field theoretical methods ||l^, |lj, |l^. Studies of integrable lattice models have added insights 



to this problem since e.g. the dependence of critical exponents on microscopic parameters and their 



behaviour due to lattice effects (back scattering, Mott transition) can be computed exactly |16, 17 



Similarly, one expects additional information from studies of lattice models for interacting electrons 



with open boundaries Q, 19, 2C]. Besides giving a deeper understanding of previous predictions 
these lattice models have features not easily included into the continuum description: local chemical 
potentials in the former lead to a sequence of bound states (see e.g. |pl[] ) which are expected to 
influence the critical properties of the boundary. 

In this paper we consider the Hubbard model on a chain of L sites subject to an additional chemical 
potential p at the first site. The Hamiltonian is given by 

L-l L , 

W = - ^ (4_,c,+i,, + /i.c.)+4n^n,tn,-i + ^iV--(iV-2iVi)-p(iVi,t+iVi,|). (1.1) 

For p = this model has been solved by Schulz using the coordinate Bethe Ansatz [Q] . Only recently, 
this solution has been extended to nonvanishing p [^] and the integrability of the model has been 
established in the framework of the quantum inverse scattering method [P2[. The Bethe Ansatz 



equations (BAE) determining the spectrum of (LI) in the A'^g-particle sector with magnetization 
M = \Ne- N^ read [|| 

g*fc,2(L+l)^^^^^.) = n^2n(%--A^)e2„(r/j+A/3) j = l,...,iVe 

/3=1 



JJe2«(AQ, -r/j)e2«(AQ, + r/j) = JJ e4tj(AQ, - A/3)e4„(AQ, + A/j) a = l,...,iV 



J = l /3=1 



(1.2) 



-ikA 



with en{x) = ^_^n , Sp{kj) = ( — ^^^7fc-~ ) ^-nd rjj = sinkj. The energy of the corresponding eigenstate 



of(0) 



IS 



£: = ^(/i---2cosfc,) +/iAr^ . (1.3) 



Using the global spin- and r/-pairing SU(2) symmetry of the Hubbard model the Bethe states extended 
by those obtained by application of the corresponding raising operators have been shown to form a 
complete basis of the Hilbert space of the system P^ j. A non-zero boundary potential destroys the 
r/-symmetry of the model and the question of completeness should be considered again. Numerical 
solutions of ( |l.2D for small p show that there exist complex combinations of two k and one A which 
coincide with one rj-pair in the limit of p ^ 0. In the following we only consider the ground state and 
the low lying excitations of the system, so we can neglect these kind of complex solutions as they belong 
to the highly excited states of the system [^]. However, for sufficiently strong attractive boundary 
potentials p > 1 we find that there exist other complex solutions which turn out to correspond to 
bound states in these potentials (note that these states do not appear in the case p < 1 studied in 
[1l9|]). These solutions need to be considered to obtain the true ground state of the system. We find 
that in spite of the presence of several complex parameters in the ground state configuration the 
low-energy spectrum of the many particle system can still be described in the Tomonaga-Luttinger 
picture equivalent to two c = 1 Conformal Field Theories. The case p > 1 will be studied in detail in 
the next section. 

2 Boundary bound states 

Prom a physical point of view it is clear that the ground state of the model contains a bound state at 
the first site for sufficiently large p. Numerical solutions of the BAE show that this is indeed the fact 
for p > 1 where a complex quasi momentum k is present in the ground state configuration. A similar 



situation has been found in the XXZ Heisenberg chain with a boundary magnetic field [21, |25[] and in 



a continuum model related to the Kondo problem |26] 



Increasing the boundary potential further we find that additional complex parameters are added 
to the gound state solution of (|1.2|). In the thermodynamic limit (L -^ cxd) we have to distinguish 
three different regions where the BAE describing the ground state are modified due to the presence 
of these complex roots:Q 



I: l<p<pi=u + VlT 



u^ 



gifc,2{L+l)^^(^^.) ^ -Q g^^(^^. _ X^)e2u{Vj + A^) , J = 1, . . . , iVe - 1 

13=1 

e2u-2t(Aa)e2u+2t(Aa) JJ e2n(Aa - r]j)e2uiK + rjj) = Y\. ^4u(Aa - A/3)e4„(AQ, + Xfs) , 



j = l /3=1 



a = l,...,N^ (2.1) 



^In principle one is free to leave the bound states empty. This gives rise to another continuum of states. These states 
become important if one considers e.g. multiple Fermi edge singularities in the presence of bound states |l3, E7|. 



with the complex solution k]\j^ = iln(p) (with exponential accuracy in the limit L — > oo) and 
t = — zsin /jjv^ = i(p — i) < li. The contribution of this bound state to the energy ( |1.3| ) is given 
by El = —p — ^ + A* — f • This complex solution corresponds to a charge bound to the first site, 
as the quasimomenta k parametrize the charge part of the states. 



II: pi <p <p2 = 2u + a/1 + 4u2 

Larger values of the boundary potential lead to an additional complex solution in the spin part: 
^N\ = i{t — u) (t > u in this region) and the following modified BAE: 

e'''^'^^^^^hp{kj) = e4u-2t{Vj)e2t{Vj) H ^2n(^i - >^i3)e2u{Vj + A/3) , j = 1, . . . , iVe - 1 

13=1 
N,~l Ni-1 

JJ e2n(Aa - ■r]j)e2u{\a + r]j) = e2t-2u{K)e6u-2t{K) Y\. '^iu{K - A/3)e4„(AQ, + A/3) , 
j=l 13=1 

Q = l,...,iV^-l (2.2) 

Again, this state can be interpreted as that of a charge bound to the surface. The physical 
excitations in the spin sector — so called spinous — correspond to holes in the distribution of 
spin rapidities A which are still real. 

Ill: p>p2 

For boundary potentials larger than the Hubbard interaction p^Au a. pair of electrons forming 
a singlet is bound to the surface, parametrized by Ajv, = siiik^^ —iu = sin/c/Ve-i +iu = i(t — u). 
The resulting BAE are 

Ni-l 

e'^^^^^+^hp{kj)e2t-4uir]j) = e2t{'nj) H ^2«(^i - A/3)e2n(r/j + A/3) j = 1, . . . , iVe - 2 

/3=1 

JJ e2„(Aa - rij)e2u{K + Vj) = Y\. ^4«(Aa - A/3)e4„(AQ, + A/3) a = 1, . . . , iV^ - 1 






(2.3) 



The energy of the second complex solution k^^-i is given by £'2 = — 2-y/l + (t — 2n)2 + /i — 2 . 

Note that region I is already realized in the ferromagnetic case with spin-| electrons only. As t = u 
(p = Pi) the index of the first factor in the A-equation of (p^) changes the sign, allowing for the 



complex A-solution. A similar change occurs in (2.1) for t = 2u (p = p2) leading to the second 
complex A;-solution. No such point exists in ( |2.3D , hence no further complex solutions are expected in 
the ground state — in perfect agreement with the physical intuition. 

Recently the BAE for the model with a boundary magnetic field {—pi{ni^ — nij)) applied at the 



first site have been constructed ||2^, ^]. This field induces an additional phase factor —e2u-2t{^a) in 



the second eq. of ( |1.2D which cancels the first factor in ( p.lj ) (up to a sign). As a consequence, we do 
not expect another complex solution to exist in the ground state besides the first one for this case. 

Using standard procedures the BAE for the ground state and low lying excitations can be rewritten 
as linear integral equations for the densities pdk) and Ps{^) of real (positive) quasi momenta kj and 
spin rapidities Aq,, respectively. Identification of positive and negative k and A allows to symmetrize 
the resulting equations with the usual result 







+ 



cos/c a2u{'n — A') 

a2M(A - i) -a4„(A - A') 



Pc 



(2.4) 



Here we have introduced ay{x) = ^ i/l+x^ ^^^ f * 9 denotes the convolution f_^dyf{x — y)g{y) 
with the boundaries ko and Aq in the charge and spin sector, respectively. The latter are fixed by the 
conditions 



^=0 _ 2[Ne-e{p-i)-e{p-p2)] + i 

CLKpc — 

■fco ^ 



Ao 



Ao 



dXps 



2[N^-e{p-pi)] + i 



(2.5) 



where 9{x) is the Heaviside step function. The driving terms of the 1/L-corrections in the different 
regions are given by 



Pc{k) = cosk a2u{ri) + 



cos k p — p^ 



vr 



7r(p^ + 1 — 2p cos k) 



e{p - pi) cos k [024(7?) + a4n-2t(7?)] 



(2.6) 



for the charge-sector^ and 



P°(A) = a4„(A) + < 



p< 1 

a2«-2t(A) +a2u+2t(A) / 

-a2t-2u(A) - a6n-2i(A) // 

/// 



(2.7) 



for the spin-sector. In terms of the dressed energies £c and Eg which satisfy the same integral equations 
as in the Hubbard model with periodic boundary conditions: 



£.9 



/i — 2 — 2 cos k 



+ 







a2u{v - A') 



(2.^ 



h I \ a2u(A — rj') cosk' —a^ui^ — A') 

the energy of the state can be expressed as: 



E 1/1 



1 



fco 



27 TT 2L 

-fco 



fco Ao 

dk e,{k)p^^{k) + f dX £.(A)pO(A) 

—fco -Ao 



+ 2L 



-(^ + - - 2) + 2e{p - i)Ei + 2he{p - Pi) + 2e{j} - ^2)^2 



+ ^'Z 



(2.9) 



Note that the index of aAu-2t changes sign at p = p2- 



3 Ground state expectation value of Ni 



The ground state expectation values for the occupation of the boundary site (A'^i) can be calculated 
from the identity (A^i) = —dE/dp. With ([2^) we obtain 



{Ni 



dk ec{k) ^J' + 



dXe,[X)^^ + 2e{p-lf-^ + 2e[p-p,f''' 



dp 



dp 



dp 



(3.1) 



-fco — Ao 

In absence of a bulk magnetic field h the ground state of the Hubbard model is known to be a singlet 
(for even particle number) corresponding to Aq = oo. In this case the system of integral equations 
(|2.8D can be reduced to a scalar one 

-ko 

with {y > 0, y + z > 0) 

e-IH 



/ko 
dk' Gll{r]-r]') COS k' e^k') 
-fcn 



(3.2) 



G^(A) = J-ReJM/f^ + f +iA 
y"- 2TTy 1 V4 4y 2y 



*'>i+'4'^' °''^^ 



2 cosh (fa;) 



(3.3) 



(^ is the digamma function). We obtain 
(Ni) = 



dE, dE2 

9{p - 1)— — e{p - p2)- 



dp 



dp 



1 



fco 



ko 



dk edk) < 



lp{k) 



p <1 



7p(fc) + I {Gl:-'\v) + Gl:+'\v)) cask I, II 



dp 

a 



(3.4) 



^ jp{k) + g- {a2tiv) - a2t~4uiv)) cos k 



III 



cos fc p +cosk—2p 



with -fpik) = ^'(^iJl^2cos^k T^ • ^^ ^^^ limit of p ^ oo only the the first two parts survive and we get 
the expected result (A^i) = 2. Some numerical results are shown in Fig. |l[ 



4 Finite size corrections 



Following [^] we can calculate the finite size spectrum of the model, reproducing the result of p9t| : 



E = Leoo + Joo + —r \ - TTT 



1 



+ 



L y 2A 2det2(Z) 

-KVs 

L 



[(Aivo - e;) z,s - (Aivf - e;) z^^f + n+ 



Ya + ^^^) [(^^° - ^^) ^- - (^^° - ^p) ^-] ' + ^^}- (4-1) 



Here A^^^^^ are non negative integers counting the number of particle hole excitations at the Fermi 



points, the Fermi velocities are given by Vc 



TTpciko) 



and Us 



4(Ao) 
TPs (Ao) 



. Z is the dressed charge matrix 



Z 



Zrr. Z, 



cc ^cs 



^sc ^ss 






(4.2) 



given in terms of the integral equation 
Cccik) Uik) \ / 1 



CcsW 6s (A) 



+ 



1 



a2«(^-A') 

a2u{^ - Vj') cos k' -a4„(A - A') 



t,cc[k ) lisc\k ) 

ecs(A') 6. (A') 



(4.3) 
The AA*"^;, are given by AA*"^ = N^ — Ln^ and AA^^ = N^ — Ln^, where n^ and n^ denote the total 



density of electrons and spin-| electrons of the reference state which we define through 



rifi 



1 fKO -I f-\Q 

J — KQ J — Aq 



(4.4) 



.(0) 



Here pcs should be computed from ( p.4D without the 1/L terms, i.e. p^^ = (note that this choice 
differs from that used in ||l^). This choice implies that for a given boundary condition AA'^^^ = 9'. 
are nonzero in the corresponding ground state. The shifts Op'^ are due to the -i--terms in ( |2.4D : 



,c,s 

p 



e; = Uj'dkp,-i + 2e{p-i) + 2e{p-p2)\ 
o; = iy_yxps-i + 2eip-pi)^ 



(4.5) 



with pc and ps denoting the solution of (2^) without the - driving term. Hence the finite size 
spectrum ( [4.1| ) determinig the bulk correlation functions |16| can be written in a manifestly particle- 
hole symmetric form by introducing AN^^ = AN^^ + 6p^, where AA'J^ denotes the change in charge 



and spin as compared to the ground state (see also [31, |3 

1 



E 



T^Vc \ 1 

Le^ + tco + — |-M + 2det2(Z) 



+ 



■KV, 



1 



1 



^KZss 



ANl'Z,, 



AN^sZcs 



AN^Z. 



+ n: 



+ n: 



(4.6) 



L { 2A ' 2det2(Z) 

These expressions simplify in certain limits (see also the corresponding discussion for the periodic 
model in ||,[l3): 

Zero magnetic field (Aq = oo): The spin part of the equations can be eliminated by Fourier trans- 
formation with the result that the matrix Z depends on the scalar dressed charge ^ = ^(fco) only [|^]: 



^cc ^cs 



^sc ^ss 



i 



which is defined as the solution of 



ko 



fco 



dk' COS k'Gllirj-rj')ak'). 



(4.7) 



(4.8) 



Furthermore, one finds the relation 9^ = ^Op, which allows to rewrite the finite size spectrum (4.1) as 



+ 



L I 24 V 2 V 



+ N 



(4.9) 



The function pc in ( [4.5D satisfies the integral equation pc{k) = pc{k) + cos k /_^ dk' G'^{'q — r]')pc{k') 
with driving term 



Pc K = - + .2,1 n M ~ ^2u[V) cos A; + cos fe ^ 

TT 7r(p^ + 1 — Ip COS /c) 



p< 1 

G^^-'\V) + Gt+'\n) I, II (4.10) 
, 024(7?) -a2t-iu{'ri) III 

The ferromagnetic case (Aq = 0): Considering the ferromagnetic case with only spin-|-electrons 

the finite size spectrum is given by 

E = Lel + fl + ^-^\-^^ + i(AiVO - ei^f + iV+} (4.11) 

and the shift O'i can be given explicitely as a function of the boundary field and the electron density 
(the Hubbard interaction is not relevant in this state): 



^'^ - = arctan ( tan — -^ 

2 TT \p-l 2 



^T,P = - o - r arctan ( ^-—^ tan ^ ) + e{p - 1) (4.12) 



5 Orthogonality exponent 

Recently, the predictions of boundary conformal field theory regarding the relation of the finite size 
corrections in the spectrum of a gapless (1 + l)-dimensional quantum system with various boundary 
conditions and scaling dimensions of certain boundary changing operators have been applied to various 
problems such as Fermi edge singularities in Luttinger liquids and the related problem of Anderson's 
orthogonality catastrophe in these systems |2|, 12, 31, ^]. Here we want to apply these ideas to study 



the second problem, namely the system size dependence of the overlap of the many-particle ground 
states corresponding to different choices of the boundary potential. For this we have to consider 
the operator Op switching on the boundary chemical potential p. Following Ref. [^] we apply the 
conformal transformation z = Le~L to get a relation of correlation functions in the infinite strip 
uj = u + iv {0 < V < L will be identified with the spatial and u with the (complex) time variable, the 
Fermi velocity is set to unity for this argument) with those on the halfplane z = t + ir, r > 0. The 
correlation function of the primary boundary operator Op in the half-plane is: 

{AA\Opin)Ol{r2)\AA) = ^^^ _^^^^,^^ . (5.1) 

Applying the conformal transformation we obtain the correlation function on the strip which is given 

by 

{AA\Op{u,)Ol{u2)\AA) ~ (-) e-^i- (5.2) 

for large Au = U2 — ui. Above we denote by \AA) the ground state of the system with vanishing 
boundary fields. The last expression can be evaluated by insertion of a complete set of eigenstates 
\BA] n) of the system with chemical potential p at the first site (boundary condition 'B') giving: 

Y,\{AA\Op\BA,n)\^e-^^S--Er\^u^^r^'^^^---^ (5,3) 



For the operator considered here the form factor (^AjOplS^d; 0) is non-zero and the exponent Xp can 
be read off to be 



L 



i,p 



- {Er - Ei"") . (5.4) 



IT 



From (^]^) we can identify Xp as the orthogonahty exponent: 

\{AA\Op\BA;0)\ = \{p\0)\r^(^^y, (5.5) 

where \p) is the ground state of the system with boundary chemical potential p. 

Using the results of the previous section we can now calculate this exponent from the finite size 
spectrum (the necessary generalization from ( |5.4D to the present case of a two component Luttinger 
liquid with different Fermi velocities in the respective sectors is completely analogeous to the one in 
the periodic Hubbard model [|l6|). The key to the correct identification of the orthogonality exponent 
is the correct choice of AA^^^ in ( [4.1[ ): as discussed above the ground state energy Eq is obtained 
by taking /S.N^s = Op=o- If we compare this energy to E^ it is crucial to compute the finite size 
corrections with respect to the same reference state. Since |0) and \p) need to be states with the same 
particle numbers N^, and Ni this implies that the correct choice of ^N^s ^^ -^o ^^ again O^Iq. 

That this choice gives indeed the desired answer is checked most easily in the ferromagnetic case: 
From ( [4.11 ) we obtain 

2 






I p=0 ~ ^^,p) — n \ ~ arctan ( tan —^ ] + — ^(P ~ 1) ) (5-6) 



approaching Xp = 2('^e — 1)^ in the limit p -^ oo. In this ferromagnetic case the many-particle wave 
function is simply a slater determinant of the one-particle functions ^k{x) oc sin(/cx) —psm{k{x — 1)). 
The product {p\0) can be evaluated numerically for finite systems leading to exponents which are in 
perfect agreement with (^). 

For the case of vanishing bulk magnetic field the finite size corrections are given by (4.9). Choosing 



AA'^^ = 0p=Q = '^dp=o = 2AA'^^ we find that there are no corrections from the spinon sector and the 
orthogonality exponent becomes 

xp = ^{e;^o-o;f (5.7) 



with 9p given in (4^). For very large p ^ oo we obtain Xp = ^^ (2 — rif^) . As we approach half filling 
Ue ^ 1 the exponent becomes Xp = -AE^^_ in Fig. ^ we present numerical data for Xp as a function 
of p for several values of rig and u = 1. 

In the general case of nonvanishing magnetic fields the exponent is given as the sum of the respective 



charge and spin part Xp = Xc + Xs with 



2^ei^(Z) [^^P=° ~ ^P^ ^'' ~ ^^P=° ~ ^P> ^''] 



X c 



^ [ {e;^, - e;) z,, - {e;^, - 0^) z,,] ' (5.8) 



2det^(Z) 



Again, this expression simplifies for p — > oo: 

,. {{2-ne)Zss- Zcsf + {Zcc- i2-ne)Zscf .^ „^ 

lim x„ = 7TZ — ^ • 5.9) 

P-oo P 2det2(Z) ^ ^ 

In Fig. y tlie exponent Xp is shown as a function of p for several magnetic fields h. 

Finally, let us remark on the effect of a second boundary potential pL at site L: the BAE (IT 



are modified by another factor Sp^{k) leading to additional shifts Op^ — Op^^Q in the expressions for 



the finite size spectrum (4.1). In this case the orthogonality exponent Xp^p^ 



-j^ \ Xp-^^pJ^ 



{piPlM ~ ( - 1 , (5.10) 

can not be obtained by simply adding the new shifts. Instead, numerical studies of the ferromagnetic 
case (see Fig. Q) suggest that the exponent Xp^p^ is given by 

^piPL ~ *^pi ' -^VL v5-llj 

i.e. the effects from the two boundaries are additive. In the framework of boundary conformal field 
theory this result is a consequence of the fact that changing the potential at both boundaries is not 
possible by the action of a single boundary changing operator Op^p^ but rather two operators Op^ 
and Op^ as becomes obvious when one switches back from the system on the strip to that on the 
half-plane (see Fig. |5|). Hence, the correlation function considered is 

\{AA\0p,{T[)0p,{n)0l^{T2)0l{T'^)\AA)\ (5.12) 

which gives (provided that |tj — r^'l <C |ti — r2|) 

\{AA\0pM)0lM)\^m{AA\0p,{n)0l^{r,)\AA)\ = ^^, _ \^^^^^^ ^^^ _ \^^^^^^ (5.13) 

for the leading asymptotic of the correlator in the semiinfinite plane. Conformal mapping of this 
expression to the strip results in (p. 11). 
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Figure 1: Ground state expectation value of A'^i as a function of the boundary potential p for u = \ 
and several electron densities n^^ (a); fixed density ng = 0.1 (b), n^^ = 0.5 (c), n^ = 0.95 (d) and several 
values of u. 
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Figure 2: Orthogonality exponent x as a function of the boundary potential p for several electron 
densities and u = 1. 
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Figure 3: Orthogonality exponent x as a function of the boundary potential p for electron density 
Ue = 0.1 and u = 1. The bold line is the exponent for the ferromagnetic case. The other ones have 
different magnetic fields h, starting with h = (dashed line) up to the critical magnetic field he 
(largest exponent for p -^ cxd) . 
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Figure 4: Numerical results for the exponent xq. 8,0.5 in the ferromagneitc case for three elec- 
tron densities and u = 1 as a function of 1/L. The plotted exponent is calculated from the ratio 

{piPL\mL/{piPLm\L+l/n.- 
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Figure 5: Conformal mapping of the infinite half plane with lu = — ln(|;) to the strip. 



